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We investigate the dynamics of a spinor Bose-Einstein condensate which is governed by an op- 
tically induced non-Abelian gauge potential. Using a ring shaped trap to confine the atoms and a 
hydrodynamic ansatz, nonlinear Josephson type equations are found to describe the system. The 
degenerate eigenstates which show rotation are solved exactly. We consider a homogenous filled ring 
and observe population dynamics between the two quasi-spin components but also space dependent 
Josephson oscillations. Stable mass currents can be observed which are induced by the constant 
non-Abelian effective magnetic field in the limit of weak interactions. For strong interactions the 
appearance of two-component dark soliton-like objects are observed. 
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I. INTRODUCTION 

The quantum ring is the Drosophila of quantum 
physics. The periodic boundary conditions have provided 
in many cases analytical solutions for non-trivial many- 
body problems and topological excitations such as per- 
sistent currents. The atomic Bose-Einstein condensate 
(BEC) [TJ [2] has in this respect proven to be a remark- 
ably versatile tool for studying a plethora of fundamental 
quantum phenomena. The Josephson effect is a prime 
example, where coherent dynamics takes place between 
two quantum states. Josephson oscillations in ultra cold 
atomic systems where indeed soon after the first realisa- 
tions of a BEC 3 , 4j [5] in the focus of interest. Character- 
ising the new condensed state lead directly to the ques- 
tion how two condensates interact. The Josephson effect 
of a weakly coupled macroscopic wave function in super- 
conductors gives a well understood archetype [TJ [6] . For 
ultracold atoms there are two different scenarios avail- 
able: a single component BEC which is separated by a 
small tunnelling barrier in a double well potential [7] , and 
a two component spinor gas trapped in a single trap with 
an internal coupling between the two condensates [5J [3] • 

The Josephson effect is a general one. ft is a mani- 
festation of phase coherence between two coupled quan- 
tum systems. Quantised vortices is another phenomenon 
which relies on many-body phase coherence, and a 
paradigm for BECs 1 . There are several techniques 
how to nucleate vortices in a BEC. With the phase im- 
printing technique, which is widely used to create dark 
solitons [in], a phase with the wanted properties is opti- 
cally induced into the gas [TTJ [I2J [13] . The most intuitive 
way, stirring with a laser beam [21 [TS] and rotating the 
trap [THl [TT1 [T5], has also been used to create vortices in 
experiments. Stirring with a laser can be a harsh tool 
if a small number of particles are used, as is often the 
case if the quantum Hall regime is the goal. Dynamic 
instabilities are the main mechanism behind the vortex 
nucleation when reaching the critical rotation frequency 
[T9j above which the vortices appear. Rotating the trap 
leads to a homogeneous effective magnetic field in the 



equation of motion of the gas in the rotating frame. A 
different approach uses optically induced effective mag- 
netic fields which is a non-rotating setup where vortices 
[101 ED HH [231 EU or solitons [25J can be created too. 
Without stirring or phase imprinting, vortex nucleation 
has also been observed in a mixing experiment of mul- 
tiple trapped condensates [25]. To date most of the ex- 
periments have been done with single component conden- 
sates, but vortices are also a feature of spinor condensates 

In this paper we investigate a spinor condensate in a 
ID ring which is subject to an optically induced non- 
Abelian magnetic field. Artificial gauge fields have at- 
tracted considerable interest lately, mainly due to the in- 
triguing analogies which they provide between ultracold 
quantum gases and magnetic effects in solid state and 
condensed matter systems. The non-Abelian gauge po- 
tentials can be created by laser assisted tunnelling [251125] 
in optical lattices, or by relying on optically induced de- 
generate dark state dynamics [30] , 

We are interested in the dynamics of a quasi-spin in 
the ring. The presence of a matrix gauge potential which 
couples to the quasi-spin provides the machinery for ex- 
otic topological states. In particular we show that the 
simplest non-trivial gauge potential which is also non- 
Abelian, and which has been used extensively in recent 
papers [3TJ [32J [33l [SI [35l [36] , can give rise to dynami- 
cally induced local spin precession. 

The paper is organised as follows. We start by intro- 
ducing the equations of motions which govern the dy- 
namics of the atomic quasi-spin, followed by an exact 
solution of the eigenvalue problem for the system. In the 
following section we use the eigenstates and a hydrody- 
namic ansatz for the BEC to derive a set of equations 
which resemble those who describe the internal Joseph- 
son effect. In Section V we investigate the spacial density 
oscillations and discuss the resulting atomic spin currents 
in the ring. Finally we conclude with a discussion about 
the influence of collisional interactions on the dynamics. 
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II. THE RING GEOMETRY AND THE ORIGIN 
OF THE GAUGE POTENTIALS 

In the following we consider a Bose gas of atoms which 
are described by a tripod level scheme [30] [37] (Fig. 
1). The tripod scheme can for instance be formed by 
three ground F = 1 states (mjr = 0, ±1) and an ex- 
cited F = tiif = state. In particular the transition 
551/2 (F = 1) <-> 5P 3/2 (F = 0) in 87 Rb or the transition 
2 3 S' 1 «-> 2 3 P in 4 He* could be used. Three lasers are 
connecting the ground-state levels tof = — 1, and 1 
to the excited state with cr + , ir, and er_ polarisations, 
with Rabi frequencies fii, 2,3(1"), and phases £1,2,3(1"). 
The resulting atom-light interaction leads to two eigen- 
states of zero energy, called dark states |£>i, 2 (r)), which 
are linear superpositions of the three lowest bare states. 
If ft = yzZi \^i\ 2 is large compared to any other en- 
ergy scale, including laser detuning, Doppler and Zee- 
man shifts and interaction energy, we can safely neglect 
transitions out of the dark state manifold. An atom pre- 
pared in such a dark state can be described by the gen- 
eral state |tf(r,t)) = ELi^M)lA(r)}, where #;(r,f) 
is the wave function for atoms in \Di). If we allow the 
atom to move, i.e., we add a kinetic energy to the light- 
matter interaction Hamiltonian, and project onto the 
dark state manifold, we obtain the effective Schrodingcr 
equation [3U] 



■ fc 5 T 



^(P-A) 2 
2m 



(1) 



where ^ T = (XP 1 , ^2) , P is the momentum op- 
erator and m the atomic mass. The key point 
here is that the spatial dependence of the laser ar- 
rangement leads to an effective matrix vector po- 
tential [35] GU SD] A„ m = iH(D n (r)\X7D m (r)). 

In addition, V m „ = (D n \Y%=i V i( T )\ D m) and 
*„ m = (h 2 /2m)(D n (r)\VB(r))(B(r)\VD m (r)) are effec- 
tive scalar potential matrices, where \B(r)) is the so- 
called bright state, i.e., the linear combination of the 
ground levels which does couple with the excited state. 
The potential Vi may include any detuning of the n-th 
laser from the resonant transition. The gauge potential 
can be shaped in quite arbitrary manner, but already a 
very simple configuration as depicted in fig. [T] leads to a 
non-trivial, non-Abelian potential [32 , 33J . 

The term p • A describes spin-orbit coupling. It will 
give rise to two single atom dispersion branches. For a 
homogeneous 2D system degenerate minima in the dis- 
persion may occur, which would for the many-body sys- 
tem cause fragmentation and preclude the formation of 
a BEC [3T]. This is not necessarily the case in a ID ring, 
as we will show later, even if the ring of course resides in 
the two dimensional plane. 

Inter-atomic interactions play an important role in the 
properties of the Bose gas at low temperatures. We will 
consider the case where the interaction energy is much 
smaller than the Rabi frequency f2, which guarantees that 




FIG. 1: (Upper left) atomic bare tripod level structure 
and (upper right) geometrical setup with three resonant laser 
beams to generate a non-Abelian gauge field in the dressed 
state basis (below). The dark state sub-manifold is encir- 
cled. The lasers drive the system with Rabi frequencies Oi = 
BaxOotter™* /\/2, H 2 = sin0 o fie iK! 7v / 2, 3 = cosB Q,e 1Ky , 

where cos 9o = V2 - 1 and Q = y / Ef=i l^il 2 - (See [32] for 
more details.) 



we always stay in the dark-state manifold. It is impor- 
tant to note that the ground states \j = —1,0,1) form 
a spin-1 Bose gas and that short-range interactions are 
dominantly s-wave. Because of symmetry collisions occur 
only in two different channels with a total spin of and 
2. These collisions are characterised by the corresponding 
s-wave scattering lengths ao and 0,2 |S]. These scatter- 
ing lengths arc in principle different, although in practice 
they are very similar. For simplicity of the discussion we 
consider ao = a 2 = a. This allows us to prevent spin 
flip collisions and we us a Gross-Pitaevskii mean field de- 
scription. The dark state spinor gas is described by a two 
component wave function 



* 2 



which gives a Hamiltonian of the form 

H=^(p-A) 2 + V(r) + ,g'(|vl/i| 2 + |* 2 | 2 ). 



(2) 



(3) 



Here, g' — 4irh 2 aN/m, where a is the s-wave scattering 
length discussed above, and N is the particle number 
with / dV(|^i| 2 + |* 2 | 2 ) = 1. The vector potential A 
can be chosen quite arbitrary. We are here interested 
in a constant, non-Abelian potential in terms of Pauli 
matrices a |3T1 [32] 



(4) 



The Pauli matrices reflect the spacial dependence while 
k is determined by the wave-vectors of the underlying 
laser fields. The Schrodinger equation for the the centre 
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of mass spinor motion then becomes with I — g'(\^fi\ 



1*2 



ih 



at 



* 2 



h 2 v 2 

2m 

h 2 K r 

m 



h A K, 2 



id x + d y ) 



h-K 

2m 
h 2 V 



(-id x - d y ) 



h 2 ^ 



2m 



(5) 




V is a constant external offset in the trapping potential 
of the two different dark states which can be controlled 
by detuning the incident laser beams from the resonant 
transitions in the tripod level structure. 



III. EIGENSTATES 



In this paper we will restrict ourselves to a ID situa- 
tion on a ring. Experimentally an effective fD ring for 
cold atoms can be achieved by strong harmonic trapping 
in the transversal direction [15] . In 2D it has been shown 
that the non-Abelian dynamics is all but trivial [3~T1 133j . 
In this work, however, we investigate a ring geometry 
where the centre of mass dynamics is ID but the corre- 
sponding spin dynamics is not, see Figure [l] There are 
a number of techniques available to experimentally cre- 
ate sophisticated ring traps [43] 04] , where even toroidal- 
shaped traps for BECs have been proposed [45] . 

We describe the ring dynamics using polar coordinates 
with the azimuthal angle ip and the radius p. The gauge 
potential in Eq. ^ is consequently transformed as 



-tiK(cr x cos ip + a y sin ip) 
-hn(—<T x sin ip + <jy cos ip) 



0. 



(6) 
(7) 
(8) 



The radial and z components are neglected in the fol- 
lowing. The ring trap provides with hu p ^ 2m(p — A) 2 
the dominating energy scale and hence any dynamics in 
these directions are suppressed by the harmonic oscillator 
ground state with a cross section s = h/mu> p . 
The relevant vector potential is now of the form 



A^ — tlHL 











(9) 
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while the momentum operator in e v -direction reads 



el t*2m t i /h 



FIG. 2: The eigenstate *+ =0 : (a) Population, (b) total (J z ), 
orbital (L z ) and spin angular momentum (S z ), (c,d) the phase 
of the spin components and (e) the gradient of the phase as 
a function of time (here for q — 0, p = 2k" 1 , k = 1). The 
white line in (e) shows the 27r-phase jump and indicates the 
quantised winding number of the second spinor component. 



pop 



(10) 



The gauge potential now depends of the angle ip. This 
(^-dependence leads to additional factors compared to the 
expression for the gauge potential in Eq. ^ . The result- 
ing Hamiltonian becomes 



9 1 

2 



V + g\^\ 



i 2 ftdtp 



2 



-if + 



<?l*| 2 
(11) 

for the energy, length in units 
of p, time in units of and g = 4 ^° jv , where N is 

the number of particles and s i = — - — \ the transverse 
cross section. The first term in the off-diagonal coupling 
vanishes in the limit p — > oo and can be seen as a kind 
of centrifugal term due to circular motion. The potential 



expressed in units of ^2 

mp 

h 



The eigenstates of the non-Abelian ring can be calcu- 
lated analytically using the Heisenberg equation of mo- 
tion for the velocity operator. The interaction term is 
merely a function of r and gives no contribution to 



v = -i[H,r] =(p-A). 



(12) 



For a homogeneously filled ring where [p, l^l 2 ] = 0, the 
Hamiltonian (111 commutes with v. The common eigen- 



states are found from the eigenvalue equation for the ve- 
locity operator which has the form 



d 

(-i- A)* = 

dip 



(13) 



with the eigenvelocity K. The Hamiltonian satisfies the 
eigenvalue equation 



E K y 



(14) 
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where the eigenenergy is given by Ek — K 2 . To solve 
the eigenvalue problem ( 13 ) we use the ansatz 



?1 



3 i(g+l)v>V|i 



</2 



(15) 



where q is an integer. Eq. (15 1 leads to a set of two 
equations for the remaining components ^ q i, 



q2> 



92 



-i/C* 9 l + (£?+ 1)* 9 2, 



with a non-trivial solution if 

K — q —ik 
in K — q - 



1 



= 0. 



The determinant gives the condition 

(K-q)(K-q-l) - n 2 = 0. 



(16) 
(17) 



(18) 



(19) 



For convenience we introduce the dimensionless quantity 
u as 



(u + l/2) 2 - l/4 = u(u+l). 



(20) 



After inserting Eq. ( |20| ) in (19 1 we get for the eigenve- 
locities 



K+ = q 



1, K_ = q — u. 



From Eqs. (16) and (21 1 we get 
1 



V2 

The eigenstates are now 



ft- 



Cv /u+ 1/2 T 1/2 e iw 
Tic^u + 1/2 ± 1/2 e^+V 



(21) 



(22) 



(23) 



where c = (u+ 1/2) -1 is a normalisation constant. From 
these expressions we see that each spinor component has 
a winding number associated to it, and that the winding 
number is always shifted by one between the two eigen- 
state components. The ratio ^ q i/^ q 2 does not depend 
on the quantum number q nor the interaction strength g 
but on the vector potential k. For any quantum number 
q the gradient of the phase difference is | VS^ — VSil = <p- 
The eigenenergies E ± of the eigenstates \E' are calcu- 
lated using (21) together with (20), 



E ± = K 2 = 



q- 



1/2 ± y/tf~ 



1/4 



(24) 



There are two branches E for every value of q which 
can be positive or negative. The states with E ± and 
q = qo are degenerate with E^ and q = — yo — 1 as 
shown in Fig. pi |\^) and |\£i 9 i) span a space of infi- 
nite many groundstates. However, the orthogonal eigen- 
states l^Sft 2 ) arc eigenstates of the total angular momen- 



form ol\ 



Hence all superpositions of the 
, /3e^|^_ x ) differ in their total angular 
momentum {a 2 — P 2 ){qo + 1/2), where the real numbers 
a, (3 obey a 2 + 1 = 1 and e 1 ^ reflects just a constant 
phase factor between the two states. 

For the ring J z is a conserved quantity as J z and H 
have common eigenstates. If the condensate is first pre- 
pared in the ring and then subsequently the lasers are 
switched on, the conservation of the total angular mo- 
mentum should prevent the system from fragmentation 




q"=4 



FIG. 3: Energy spectrum of the ring. The dashed (red) 
curves belong to q + > and are degenerate with q~ < 0. 
The solid (blue) curves belong to q~ > and q + < 0. The q 
is the states winding number and the ± indicates the branches 
{E±. All energies are two times degenerate. Superpositions 
sin^M^,) + cos Sj4'i 90 _ 1 } have different J z as shown in the 
insert for go = 0. 



IV. THE INTERNAL JOSEPHSON EFFECT 

In the following we will use a hydrodynamic descrip- 
tion of the condensate [UJ to obtain a set of equations 
describing the dynamics on the ring. We write an ansatz 
for the wavefunction as 



*J2 



,iSi 



-Vn 2 e 



iS 2 



(25) 



with the (normalised) local density rii(ip) and the phase 
Si(<p). 



By using the transformation in Eq. ( 25 1 together 



with the Hamiltonian from (111, the time dependent 



Schrddinger equation is rewritten by multiplying with the 
hermitian conjugate Vtj from the left, 



- iTft 



(26) 



turn operator J z 



S' 



-icL 



From Eq. ( 26 ) we get with a straightforward calculation 



ha z as they fulfil (see appendix \m after separating in real and imaginary 
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<b> t * 2m p7 h 



u 



<e) C2mn 2 /h 
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(fl t*2m(j z /h 




<g) t * 2m || 1 h 



FIG. 4: The Internal Josephson effect by applying an external 
potential V = 0.5: (a) The population of the two quasi-spin 
components is shown. The black and solid curve shows the 
population simulated numerically while the red and dashed 
line is the analytic approach from Eqs. (40 1, (b) shows total 



(J z ), orbital (L z ) and spin angular momentum (S z ). In (c,d) 
we plotted the density of the BEC, the minima are black 
and the density maxima are white. (e,f) show the phase and 
(g) the phase gradient VS2 of the second component for the 
AC- Josephson effect. In all cases p = 2k~ with an initial 
eigenstate (c^/u exp(0), —icy/u + 1 exp(i<^) ) 



parts the coupled nonlinear equations of motion for the 
phases Si(ip,t) and the local densities rii((p,t): 



Si = [ 



1 



-^r ^ " '2 ( ^ 5i) 

-(K 2 + l/ + . 9 (ni+n 2 )] 

mi \ 2 2n 2 



msA + ^sinA)(27) 



u x = 



( - d v S\d v nx - rii(9^Si)^ + 2K % /n 1 ?i 2 

x((d v S 2 - ^)sinA - ^^-cosA 
V 2 2n 2 



(28) 



So 



-(d v n 2 ) 2 , 9>2 (d v S 2 ) 5 



^_ 

8n| ^ 2n 2 2 

-(K 2 -y + . 9 (n 1 +n 2 ))] 

+ a„S 1 )co8 A sin A)(2<>) 

^/n 2 V 2 2rii ' 



h 2 = ( - d v S 2 d v n 2 - n 2 (d 2 S 2 )j + 2ny/nin 2 

*({d v Si - ibsinA- %^1 C osaV (30) 
V 2 2ni / 

where we have introduced A = S 2 — ip — Si . For a com- 
pletely filled ring with a constant density, d v rii = 0, the 
equations simplify considerably, 



Si 

fix 
S 2 

fl2 



1 
2 



(9^5i) 2 - (k 2 + V + g(m + na)) 



(cLS 2 - -) cos A 

Hi 2 



(31) 



-ni(a;5 x ) + 2^7711^2 (-- + <9 V S* 2 ) sin A(32) 



■ 2 (5,S 2 ) 2 -( K 2 -K + 5 (n 1 + n2 ))] 



^4 + Vu)cosA 
n 2 2 



(33) 



-n 2 (d;S2) + 2 K% M^(-- + fl^Si) sin A(34) 



The cigenstates (23) can be seen to fulfill the relations 



S2 — <p — Si =0 and d 2 Si = 0, hence cos A = 1 and 
sin A = 0. The time evolution of the densities is then 
given by 



1 



ni 



-2«v/nin2(- - d v S 2 ) sin A = 0, 



h 2 = -2K % /n i n 2 (- - d v Si) sin A = 0. (35) 

The phase difference A = S 2 — ip — Si, on the other hand, 
evolves according to 

Si-S 2 = l -{{d v S 2 ) 2 ~ {d v S!) 2 ) -IV (36) 



n 2 , 1 



d v S 2 ) 



(± + d v S 1 )-(± + d v S 1 ) = 0, 



(37) 



where we have used the eigenstate properties d v S 2 — 

l + ^5iand^-^=iandV = 0. 

Oscillations which are driven by the phase difference 
between two quantum fluids where first investigated in 
solid state physics using superconductors with a small 
tunnelling barrier, e.g., a grain boundary. The predic- 
tion of AC and DC currents in so called superconducting 
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FIG. 5: Internal Josephson effect for non eigenstate 
(0.5,-0.5iexp(ivp) ) T : (a) The quasi-spin population 
dip |*i| 2 , i = 1,2; (b) the total (J z ), orbital (L z ) and 
spin angular momentum (S z ) are plotted for each quasi-spin 
component, (c) and (d) show the density of component 1 and 
2, where the minima are black and the maxima in white, (e) 
and (f) is the the phase of component 1 and 2, ranging over 

±7T. 



which becomes 

hx(tp) = K^/n\ri2 sin[wi], 

"■2(v) = —K^Jri\ni sax\iot\ = —fix. (39) 

These coupled differential equations can be solved exactly 
with =p for n 12 respectively, 



nx,2(t) = (y "1,2(0) cos[k/o; sin 2 [cji/2]] 



T ^2,1(0) sin^/w sin 2 [c^/2]]) . (40) 



Eq. (40) describes population oscillations with a fre- 
quency of lo = 28V. The observed frequency in Fig [4] 
is only slightly smaller then lo = 25V frequency. 



V. THE SPACIAL JOSEPHSON EFFECT AND 
VORTEX NUCLEATION 

In solid state materials the phase difference between 
the two macroscopic wave functions is denned as the dif- 
ference at the boundary surface of bulk materials. Each 
phase is considered constant near the boundary surface 
and the phase inside the bulk material is assumed not 
to contribute. The situation in the ring offers a different 
physical scenario as the relevant phase difference is not 
necessarily the same in every point of the ring. 



interference devices (SQUIDS [IS]) was a great success of 
the BCS theory of cooper pairs. SQUIDS are nowadays 
widely used e.g. in magnetic field sensors. The set of 



equations ( 35 ) and ( 36 1 describing the non- Abelian ring 



resembles the Josephson equation and does indeed give 
rise to coherent effects like population oscillations in the 
ring. 

If we apply a constant potential ±U to the two dark 

q2 



states, or if we change the population ratio 



of an eigenstate (22 1, the system shows population os- 
cillations as illustrated in Figs. [3] and [5] for a non- 
interacting system with g = 0. The total angular momen- 
tum J z is always conserved, whereas the oscillations of 
L zi {t) = (9i(t)\L z \9i(t)) and (9i(t)\S a \9 t {t)) for i = 1,2 
stem from the population oscillation of (^i(t)\^i(t)). 

The population oscillations' frequency can be esti- 
mated for a small perturbation V = SV of an eigenstate 
^I /± as the phase difference Si — £2 in Eq. (36 1 is no 
longer constant, 




Sx - S 2 = -2 SV = 



(38) 



FIG. 6: Spacial Josephson oscillations: Single component 
density (a, b) and overall density (c) showing spacial oscil- 
lations for the initial state (1,— i). The spin orbital angular 
momentum S z (ip,t) oscillates locally, but not globally. (Com- 
pare with S z (t) in Fig. [7](f)). 

To illustrate this we consider an initial state of the 
form 



With A = Lot we get for the lowest eigenstate (q = 0) 
Josephson population oscillations as seen from Eq. ( 35 ) 



^ = 



1 / 1 



4tt \ - 1 



(41) 
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which leads, for a non-interacting gas (g 
ics according to (31 1-( 34 ) 



2 «V"2 
-K + COS ip, 

hi(f) = —Ky/niU2 sin tp, 



Si 



■ 2 rc-v/ni 



) = —K\Jri\n-2. sin y>. 



0), to dynam- 

(42) 
(43) 
(44) 
(45) 



These equations have now a nontrivial spacial depen- 
dence where the density changes with same sign for both 
components of the system. As we have rij ~ sin (p we 
expect a spacial distribution with one maximum and one 
minimum on the ring. Moreover, the time derivative of 
the velocity field is d v Si = — cLSa if n i — n 2- Fig. ^ 
and [7] show the evolution of (41 1 in time. Initially we 
see the dynamics described by Eqs. (42|-(45| with spa- 
cial sinusoidal density oscillations rij ~ sin <p and currents 
d^n^ip) = -d<f,n 2 (<p). 

For longer times, the ring is described by the coupled 
and nonlinear eqs. (27l-(34l and we observe not only the 



spacial density distributions but also counter propagating 
currents. For the chosen gauge potential the effective 
magnetic field stems from the non-Abelian part and is 
given by a constant field in the z-direction 



B = V x A - —(A x A) = hn 2 
2/T ; 





-1 



(46) 



where both spinor components start to rotate against 
each other. 

In the phase-time diagram for S(tp,t) in Fig. [7] the 
moment of establishing (and vanishing) a quantised cir- 
culation in each spin component can be found where 
the phase singularities appear. In the velocity field 
Vi = d v Si(ip,i) we can clearly observe regions with posi- 
tive and negative velocities v^. 

Persistent currents and vortex nucleation have been 
widely discussed in literature [HI [TBI P^l W?\ TO P^l 
l20l I26} |27| 149] as the nucleation of rotation needs typ- 
ically a dynamic instability. The two component spinor 
gas offers the option, that the components rotate against 
each other. As the gas' topology has to be preserved, 
the rotation emerges in counter rotating pairs. The nec- 
essary break up of the coherent wave function happens 
at the single component density minima, ^(tp)! 2 = 0, 
where it costs no energy to change the phase. Never- 
theless, the condensate itself does not break into parts 
as the singularities in the phase-time diagram appear 
simultaneously but spatially separated (see Fig. [7] (e)) 
by approximately 15°. In addition, the orbital angular 
momentum L z i — (^j(t)| — i9 v |^i(t)) changes smoothly 
for each component in time and has no discontinuities. 
It is still possible to define a preserved winding num- 
ber q by counting the ±2tt phase jumps with the proper 
sign for both spinor components. The spin orbital mo- 
mentum, S z (t), is conserved in time, but not locally as 
|*i(y>,t)| 2 - \^ 2 {(p,t)\ 2 ^ 0, which is shown in fig. |6](d). 
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FIG. 7: The phase and the phase gradient of each spinor 
component are depicted as function of time. The singularities 
in the phase-time diagram (a),(b) indicate vortex creation and 
annihilation. The gradient of the phase is overall positive 
(negative) sign once the current is established (white / black 
area in (c),(d)). The black and white lines are the ±2ir phase 
jumps, (f) The single component orbital angular momentum 
L z i(t) oscillates while the total spin angular momentum J z (i) 
is constant. 



VI. THE ROLE OF INTERACTIONS 



As shown in Section |ITT| the eigenstates do not depend 
explicitly on the interaction strength g and also the inter- 
nal Josephson oscillations are not qualitatively affected 
by the interactions. In contrast the spacial oscillations 
and thereby especially the mass currents are much more 
fragile in the presence of interactions. This is indeed to 
be expected. The elementary excitations of a standard 
single component condensate generally also depend on 
the interaction, which for instance give the expression 
for the speed of sound in the gas. An additional restric- 
tion which needs to be kept in mind for the spin-orbit 
coupled gas is the size of the chemical potential, which, 
in order to fulfill the adiabatic approximation so that the 
atoms remain in the dark state manifold, must be smaller 
than the total Rabi frequency. 

It is first of all important to note that an arbitrary 
initial state can induce large amplitude oscillations and 
population transfer between the spin components. In the 
presence of interactions, and consequently nonlinearitics, 
the dynamics is inevitably going to be complicated if the 
initial state is not similar to one of the eigenstates. This 
is indeed always the case, also for single component con- 
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densates, however, the spinor character of the ring stud- 
ied here tends to enhance this effect due to the increased 
parameter space. 

fn order to prepare a well defined coherent spinor con- 
densate as initial state several criteria needs to be ful- 
filled. First of all, as mentioned above, the adiabatic- 
ity needs to be preserved. Secondly, the spin-orbit cou- 
pled gas is notoriously problematic due to fragmentation 
which would preclude the formation of a coherent con- 
densate |31j . To circumvent this problem one can pre- 
pare a well defined condensate in one of the three lowest 
to = —1, 0, +1 states, and then transfer adiabatically the 
condensate to the dark state manifold. By doing so the 
coherence is expected to be preserved, and a well defined 
Gross-Pit aevskii equation can be used to describe the 
gas. 

The resulting dynamics would however typically be 
well within the nonlinear regime, meaning that a 
Bogoliubov-de Gennes treatment, or in other words, a 
linearisation of the Gross Pitaevskii dynamics, is not ap- 
propriate. This is because the adiabatically transferred 
state is not necessarily going to be a steady state solu- 
tion in the new spin-orbit coupled setting. Any attempt 
to identify eigenmodes and eigcnfrequcncics will therefore 
inevitably be strongly influenced by the collisional inter- 
actions. In order to analyse the elementary excitations 
using the Bogoliubov-de Gennes treatment one should be 
close to the ground state of the system and study small 
deviations from it. But this is not necessarily straight- 
forward because the fragmentation will prevent us from 
preparing, i.e., cooling, the condensate directly into a co- 
herent ground state. Therefore, in most cases, there will 
be dynamics taking place, but some adiabatically pre- 
pared initial states will be better approximations of a 
coherent steady state than others. 

We have simulated the dynamics of the ring conden- 
sate using a single scattering length g. Strictly speaking 
scattering takes place between the bare electronic levels 
and not in the dressed state basis of the dark states 
A more detailed analysis has to take different scattering 
lengths between the atomic hyperfine levels of the tripod 
system into account and project these on the dark state 
basis. However, if the scattering lengths are close to each 
other, one can justify the usage of a single parameter g 
[50] . For 87 Rb the scattering lengths are typically within 
a few percent |51j . 

Our simulations show in Fig. [8] a significant influ- 
ence on the mass current creation above an interaction 
strength for which g > 5. Instead of inducing counter 
propagating currents, the spacial Josephson oscillations 
lead to several pairs of soliton like objects which circu- 
late around the ring [52] . The non- Abelian magnetic field 
is not strong enough to induce stable currents any more 
(see Fig. [9). 
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FIG. 8: Strong Interactions destroy the mass current. (a,b) 
Total, orbital and spin angular momentum, (c,d) density, and 
(e,f) phase gradient with interactions for g = 2 (left) and 
g = 25 (right); Different than in the interactionless case (g = 
0) a regime with persistent currents over the full ring is not 
established. The initial state is l/v / 47r(l, — i) T . 



VII. SUMMARY 

The bosonic non- Abelian ring offers a rich toolbox for 
cold atom physics. We investigated a homogeneously 
filled ring which can be fully described and understood 
in terms of the well known internal, nonlinear Joseph- 
son effect. Moreover, it was found that this mechanism 
can cause large density oscillations and give rise to grey 
soliton-like structures. 

The eigenstates were calculated exactly. They show 
mass currents where detunings in the external trapping or 
deviations from the eigenstate population ratio leads to 
the AC- Josephson effect. Spacial Josephson oscillations 
appear for the most simple configuration of an equally 
filled ring without any initial rotation. Already the non- 
interacting case shows non-trivial behaviour as spacial 
density oscillations and oscillating mass currents appear. 
In the interacting case the coherent oscillations in the 
quasi-spin populations and the spin currents breaks down 
for sufficiently large interactions, which results in the ap- 
pearance of soliton-like structures. 

Observing the spacial density oscillations would exper- 
imentally validate the non-Abelian nature of the opti- 
cally induced vector potential. In the dilute regime the 
nucleated rotation vanishes periodically but the found 
effect might be used to initiate rotation in the system 
and switch off the laser field in a coherent way to pre- 
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FIG. 9: Density dynamics for g — 25 (a,c) and g = 150 
(b,d,e,f). Phonons can be seen to propagate at the speed of 
sound (middle), while the grey soliton like objects propagate 
considerably slower. In the strongly interacting regime the 
two condensate components start to separate (e, f). 



serve the rotating states in their traps. Furthermore we 
found a mechanism to nucleate persistent currents with- 
out any dynamical instability or atom-atom interaction 
necessary. 

Finally it is worth noting that there is extensive dis- 
cussions in the literature about non-Abelian Josephson 
effects, which is a field theoretical generalisation of the 
conventional Josephson effect and can be connected to 
the Higgs mechanism [53 54 . In our paper, we use a non- 
Abelian vector potential to couple the two wave function 
dynamically, but the observed spacial Josephson effect is 
still a realisation of the Abelian Josephson effect. For the 
non-Abelian effect the setting is, e.g., a double well po- 
tential separating an Abelian and a non-Abelian region. 
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with the hydrodynamic description of the condensate 
wave function, i , ^2 ) T = (nie lSl , — in2e lS2 ) T . For the 
left hand side of (Al) we use the relation 



= i ^^ n i - V^^^A (A2) 
ot 2^/nlot at 

After multiplication with (^\, from the left the time 
dependence reads 



#*id t #i = i l -ni - n l S l 



(A3) 



On the right hand side of (All we have diagonal and 
off-diagonal parts. The diagonal parts can be calculated 
with 



8n 4 



(d v n i y--(d'n i ) + ^(d 9 S i y. (A4) 



For the off-diagonal terms, mixing the two components 
^2, we use 

y/^e- iS *d vy fce iSi = (A5) 
e i( Si -s s ) 1/2 [™I d(pni + i^^ j e i( - s '- s ^d v S l . 



The first off diagonal matrix element has the form 



r — 1 

T 

1 i 

2 ~ 2n 2 



d v ii 2 + d v S 2 



(A6) 



while the second off diagonal element reads 



APPENDIX A: DERIVATION OF JOSEPHSON 
EQUATION 



1 

2 

1 i 

2 ~ 2r^ 



d v ni + d v S\ 



(A7) 



The set of nonlinear coupled equations ( 27 1 to ( 30 1 is 
obtained from ( 26 1 



Collecting the real and imaginary parts of (A3 1, (A4| 
( A6 1 and ( A7 1 gives using e 1 * 



i sin (f>+cos $ the coupled 
inear equations for the phases Si(ip,t) and the local 



(Al) densities n,i(<p, t). 



